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. . ABSTRACT 
<_ 

I Using deep J, H and JsTs-band observations, we have studied the near-infrared 

extinction of the Nuclear Bulge and we find significant, complex variations on small 
physical scales. We have applied a new variable near-infrared colour excess method, 
V-NICE, to measure the extinction; this method allows for variation in both the 
Q ■ extinction law parameter a and the degree of absolute extinction on very small physical 

' scales. We see significant variation in both these parameters on scales of 5". In our 

C/3 , observed fields, representing a random sample of sight lines to the Nuclear Bulge, we 

measure a to be 2.64 ± 0.52, compared to the canonical "universal" value of 2. Our 
measured levels of are similar to previously measured results (1 ^ ^ 4.5); 
however, the steeper extinction law results in higher values for Aj (4.5 ^ Aj ^ 10) 
^ ' and Ah (1-5 ^ Ah ^ 6.5). Only when the extinction law is allowed to vary on the 

t — I smallest scales can we recover self-consistent measures of the absolute extinction at 

OO . each wavelength, allowing accurate reddening corrections for field star photometry 

0^ ' in the Nuclear Bulge. The steeper extinction law slope also suggests that previous 

conversions of near-infrared extinction to Ay may need to be reconsidered. Finally, 
we find that the measured values of extinction are significantly dependent on the filter 
' transmission functions of the instrument used to obtain the data. This effect must be 

0^ . taken into account when combining or comparing data from different instruments. 
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1 INTRODUCTION be 4v ~ 30mag (Rieke & Lebofsky, 1985; Catchpole et al., 

^ , , , , „ n , >T , 1990;Blumetal., 1996; Dutraetal., 2002, 2003, and others). 

Study or the structure and stellar content of the Nuclear mi i- .li • i- • j_i ^ i mnv r ^^ 

.. r , 1 he reason tor this discrepancy is that only ~ 1070 ot the m- 

Bulee IS extremely diihcult as it is one ot the most highly , , ,i i. . i. . . i i i i . 

,,vTn -I-,! terstellar medium is distributed homogeneously throughout 
obscured regions or the Galaxy. The Nuclear Bulge is es- at i n i i • • r,r,iw • j_ • i 

° ° the Nuclear Bulge and the remaining ~ 90% is contained 

timated to contain a narrow layer of molecular hydrogen -^i • i i i ^ i j_ i-^i at i 

r'. '^ ° withm dense clouds falling only a tew percent ot the Nuclear 



Bulge (Catchpole et al., 1990). 



throughout its extent with scale height /i ~ 30 — 50 pc and a 
total mass of M ~ 10* Mq (Mezger et al., 1996; Rodriguez- 
Fernandez & Martm-Pintado, 2005). However, this accounts 

for only half of the extinction between the observer and the , „, t 

1.1 The Extinction Law 

Galactic Centre as it is estimated that there is an equal 

amount of material between the Sun and the Nuclear Bulge The relationship between extinction and wavelength is gov- 
as there is within the Nuclear Bulge in the direction of the erned by the composition of the material causing the ex- 
Galactic Centre (Gordon et al., 2003). tinction; in most cases, this is a diffuse mixture of gas 
If this material were distributed homogeneously over and dust (small grains and molecules; Fitzpatrick, 2004), 
the sky, it would result in visual extinction of J5 100 mag. and is known as the "extinction law". In the visual wave- 
Numerous authors have measured the extinction in the vi- length regime, the extinction law is generally expressed as 
sual and near-infrared and found the level of extinction to the ratio of absolute to relative extinction Ry = Ay / Eb-v ■ 

In the near-infrared wavelength range the extinction law 
is described by a power-law relationship, A\ oc A~", and 
* e-mail: andrew.gosling@oulu.fi can be expressed either as a ratio of relative extinctions 
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Ej-h/Eh-k or, as the slope of the power-law a (Rieke 
& Lebofsky, 1985; CardelU et al., 1989; Mathis, 1990, and 
others) . 

The original, and generally applied, extinction law for 
the diffuse interstellar medium (i.e.: not dense, known molec- 
ular clouds) was calculated in the 1980s using relatively 
primitive (by today's standards) infrared detectors and us- 
ing a small sample of individually selected stars whose spec- 
tral types were well known (Savage & Mathis, 1979; Rieke & 
Lebofsky, 1985; CardeUi et al., 1989; Mathis, 1990; Catcli- 
pole et al., 1990; Whittet et al., 1993). Based on the small 
sample of sources for which they had spectral observations, 
they expressed the ncar-infrarod extinction law as "univer- 
sal" , a notion that has been adopted in almost all studies 
since. There have been many values given for the extinc- 
tion law, with values of Ej-h/Eh-k ranging from 1.28 to 
2.09 (Davis et al., 1986; Jones & Hyland, 1980); however, 
this variation has been attributed to the different magni- 
tude systems being used (see Kenyon et al., 1998, for an 
example of photometric systems affecting the extinction law 
determination) . Up to this point the possibility that the un- 
derlying extinction law could vary on small scales had not 
been considered extensively. 

However, in the last decade, as infrared detector tech- 
nology has advanced, it has been possible to obtain much 
higher-resolution, deeper surveys of the sky in the infrared 
band {2MASS, DENIS, UKIDSS, VISTA). Based on these 
observations, it has become evident that the extinction law 
is not "universal" as previously thought, but is in fact highly 
variable from point-to-point (Kenyon et al., 1998; Udalski, 
2003; Messineo et al., 2005; Nishiyama et al., 2006; Fitz- 
patrik & Massa, 2007; Froebrich et al., 2007). This varia- 
tion is most evident towards the Bulge and in denser re- 
gions such as molecular clouds (Draine, 2003). Nishiyama 
et al. (2006) have shown there is variation in the extinc- 
tion law within their observations of the Nuclear Bulge. 
They measure Ej-h/Eh-k = 1-72 ± 0.04, equivalent to 
an a = 1.99 ± 0.02 for the outer regions of the Bulge, but 
they find a spatial variation with values of a of 1.96, 1.97, 
2.09 and 1.91 for the North-East, South-East, North-West 
and South- West regions respectively. Froebrich & del Burgo 
(2006) demonstrated how the assumption of a power law 
extinction-wavelength relationship can be manipulated to 
provide equations to measure that relationship given pho- 
tometry in three wavebands; this is the method we will uti- 
lize in this paper. They applied this method to the Galactic 
anti-centre using data from 2MASS (Froebrich et al., 2007) 
and found a large spread in the measured extinction law pa- 
rameter a. However, when calculating the absolute extinc- 
tion, they took an average of the spread in extinction law 
parameter to convert from colour excess to extinction. Most 
recently, Nishiyama et al. (2008) showed in observations of 
the outer regions of the Galactic Bulge that the conversion 
factor between near-infrared and visual extinction was much 
higher (Av/Ak ~16, compared to the canonical value of 
~9) than previous studies implied, indicative of a steeper 
extinction law slope (a higher value of a). 

1.2 Previous Extinction Surveys 

The first large scale extinction mapping of the Nuclear Bulge 
in the near-infrared was undertaken by Gatchpole et al. 



(1990) using H ~ K and if mag from observations on the 
1.9 m telescope at SAAO Sutherland. They measured the 
modal colours and magnitudes from histograms of all the 
stars in in 400" x 400" cells, and compared these values to 
the giant branch of 47 Tuc to measure 5 ^ Av ^ 35. They 
observed a general trend of increasing extinction levels to- 
wards zero latitude in the plane, and also made a "hand- 
drawn" map of the distribution of dark clouds obscuring all 
background sources in their observations. Additionally, in a 
number of the histograms of colours and magnitudes of the 
stars, they observed a double peak in the distributions indi- 
cating more than one source of extinction along the line of 
sight. 

Schultheis et al. (1999) used data from the DENIS sur- 
vey of the Bulge to form an extinction map with a resolu- 
tion of 1'. They used a two-to-one over-sampling, so for their 
resolution they were sampling stars within a 2' radius. They 
used the extinction law of Glass (1999) to convert their mea- 
sured colours to absolute extinction. They measured Ay as 
high as 30, although in many regions this was only a lower 
limit as in high extinction regions the majority of sources 
were obscured, especially in J. They also reported the pres- 
ence of a clumpy and filamentary distribution to the areas 
of high extinction, especially close to the Plane. 

The last large-scale extinction mapping of the Bulge 
before this work was performed by Dutra et al. (2003) over 
10 dcg X 10 deg using data from 2MASS. They sampled 4' x4' 
regions and measured extinction ranging from Akq = 0.005 
at edges of the region to Ak^ = 3.2 close to Galactic Centre. 
As with the Gatchpole et al. (1990) map, they observed a 
general reduction in the level of the extinction away from 
the plane, and also a filamentary structure to the regions of 
highest extinction. 



2 EXTINCTION THEORY 

As discussed in Section 1, the extinction in the near-infrared 
is generally accepted to have a power law relationship to 
wavelength: 

Aa = C7xA-". (1) 

where C is a dimensionless constant. We use a similar 
method for measuring the extinction law and the value of 
absolute extinction as outlined and discussed in Froebrich & 
del Burgo (2006). Colour excess, {E{\\ — A2)), for an object 
is the difference between the observed colour for that object, 
and its intrinsic colour if there were no extinction: 

{£(Ai - A2)) = iyrnf"' - mf"'') - (mi°* - m2°*) . 

As the observed magnitudes are a function of the in- 
trinsic magnitude and the extinction at that wavelength, 
rn°^^ — m'"' + A\, the colour excess becomes 

(£(Ai - A2)) = ^Ai - Ax^ 

which allows the colour excess to be expressed £is a function 
of two wavelengths and the dimensionless constant C 

{S(Ai - A2)> = C x (Ap - AD (2) 

Using Eq. 2 for two different colour excesses with a com- 
mon wavelength, it is then easy to express the ratio of the 
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two different colour excesses as a function of only wavelength 
and the extinction law slope a. 



{E{\,-\2)) _ {^y 



{E{\2-\z)) 



(3) 



By comparing the measured near-infrared colour excess 
variation across images, the extinction law can then be cal- 
culated at each spatial position, the value of C can then 
be derived for each location using Eq. 2 and thus the ab- 
solute value of extinction for a certain wavelength can be 
calculated using Eq. 1. This is a more detailed (in that it 
considers the extinction law a variation as well as absolute 
extinction) and variable version of the NICE and NICER 
extinction calculation methods of Lombardi & Alves (2001) 
which we term V-NICE. 



3 DATA 

This work uses data from deep imaging of the Nuclear Bulge 
of our Galaxy, within a region 1.6° x 0.8° in (/, h) centred on 
Sgr A*. We observed 26 locations within the Nuclear Bulge 
using the ISAAC camera on the VLT (Bandyopadhyay et al., 
2005). The 26 fields are located in areas away from known 
regions of star formation or other known structures within 
the Nuclear Bulge. The ISAAC instrument has a 2.5 x 2.5 
arcmin^ field of view with a pixel resolution of 0.1484". 
The observations were obtained on nights with <0.6" seeing. 
Each field was observed for a total of 360 s in each filter, J, 
H and Ks- A random 20" offset was used between integra- 
tions using the standard ISAAC jitter templates to allow 
self-flattening of the images. The total 6 minutes integra- 
tion per filter gives limiting magnitudes of J=23 (S/N=5), 
H=21, and 7^s=20 (S/N=10). 

Initial reduction (flat-fielding, removal of bad pixels, 
and sky subtraction) was performed with the ESO/ISAAC 
pipeline reduction software. These image products were then 
astrometrically locked to the Two Micron All-Sky Survey 
{2MASS). Source positions and magnitudes were extracted 
using SExtractor (version 2.3.2; see Bandyopadhyay et al., 
2005, for further details). We used sources that were com- 
mon to both the SMASS catalogue and our own observations 
to match the photometry of our sources to that of 2MASS. 
This linear conversion to the 2MASS photometric system 
corrected for slight variations in the photometric calibration 
supplied by the ESO/ISAAC pipeline and allows comparison 
of this work to the previous extinction mapping of this region 
which was based on 2MASS data (Dutra et al., 2003; Froe- 
brich et al., 2007). The photometric and astrometric data 
for the sources were combined into a series of catalogues, 
one for each of our 26 VLT flelds. 



4 CALCULATING THE EXTINCTION 

We generated a grid of positions across the flelds at each of 
which the extinction is calculated. The grid has a 5" sep- 
aration between points and the grid is defined around the 
centre of each field. At each grid point, all stars within 10" 
and 20" are extracted from the source catalogues and the 
median value of the stars' colours (magnitude of the star in 
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Figure 1. Histograms of the value of the extinction law param- 
eter a, calculated from the measured colour excess ratios using 
Eq. 3 (red). The best fit single gaussian to the distribution has 
H = 2.64 with a = 0.52 is overlaid on the data (black). 



a particular filter e.g.: J) and colour- index (difference be- 
tween the brightness of a star in two colours e.g.: H — Ks) 
were calculated for that grid point. The median is used as it 
reduces the influence of foreground stars on the values cal- 
culated (the effect of the use of the median instead of the 
mean is discussed in Froebrich et al., 2007). 

Median colour excess is then obtained by subtracting a 
value for the average intrinsic colours for the stars. We used 
the average of the colours of G, K and M giant stars as taken 
from Cox (2000, and references therein) to obtain colour 
excesses for the Nuclear Bulge stars. We use G, K and M 
type stars as these are the most common stellar types in the 
Nuclear Bulge. We also set a threshold limit of 5 stars within 
the 10" sampling box as a minimum for the number of stars 
required for the scripts to calculate a speciflc extinction. If 
there are fewer than this number in the sample region, the 
value for that point for both the 10" and 20" sampling is 
set to the average value measured for the whole fleld. 

The colour excess ratios for each grid point were then 
used to identify the value of the extinction law parameter 
a by comparison to tables of values generated using equa- 
tion 3. We then used this value of a to calculate the value 
of C at that position using equation 2, and then the abso- 
lute extinction in J, H and /fs-bands was calculated using 
equation 1. 



5 THE EXTINCTION 

5.1 Measured Extinction Law 

Figure 1 shows a histogram of the values of a measured for 
all the flelds using the colour excess ratio method in Eq. 3 
with a single gaussian flt to the distribution. 

The distribution of a is well flt by a single gaussian with 
^ = 2.64 and a = 0.52. Errors in the intrinsic measured 
magnitudes of the stars are the largest source of error in the 
measurement of the colour excesses from which each value 
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Figure 2. Histograms of the values of absolute extinction calculated for all grid points. Each histogram represents the absolute extinction 
of a specific wavelength; left is Aj, centre is Ajj and right is A^g. 



of a is obtained, but the large sample (> 30, 000 positions) 
results in a vanishingly small error in the mean and standard 
deviation as measured from the histograms. The binning for 
the histograms is 0.05 in values of a; the error in the fit of the 
gaussian to the histogram is comparable to one histogram 
bin. 



5.2 Measured Extinction 

Using the value of a measured at each point in the fields, we 
calculated the absolute extinction from the colour excesses 
using equation 2 to obtain a value for C and equation 1 
to calculate the extinction for each wavelength. Figure 2 
shows histograms of the values of absolute extinction for 
each wavelength, calculated at all grid points in the 26 VLT 
fields. 

All three distributions of absolute extinction have com- 
mon features. The values of extinction are generally confined 
to a single range of values, with almost no cells having val- 
ues outside this range. Within this range, there is a small, 
low-extinction component, the main component of the ex- 
tinction is then one large peak at an absolute extinction a 
little higher than the low extinction component. Beyond this 
initial peak there is a plateau extending to higher levels of 
extinction before a relatively sharp cut-off. At the cut-off 
there is then a smaller component that forms a tail extend- 
ing to very high levels of extinction. That these features are 
common to all three wavelengths suggests that the differ- 
ent wavelengths are showing different representations of a 
common distribution of extincting material and not differ- 
ent materials acting at different wavelengths. 

The majority of the Aj extinction is in the range 4.5 ^ 
A,j ^ 9, consisting of three peaks. The low-extinction com- 
ponent of the distribution peaks at Aj = 4.8, the main com- 
ponent of the distribution peaks around Aj = 6.1 and the 
higher extinction plateau region is centred around Aj — 7.8. 

The Ah extinction distribution is in the range 1.5 ^ 
Ah ^ 6.5 with a proportionally smaller low extinction com- 
ponent and larger high extinction component. The low ex- 
tinction component peaks around Ah ~ 2.1. The stronger, 
main peak is around Ah = 2.5 and the higher extinction 
plateau region stretches over about 1.5 mag in Ah centred 
around Ah = 3.5. There is also an small, even higher ex- 
tinction tail extending between 4.5 ^ Ah ^ 6. 




4 6 8 

A (magnitudes) 

Figure 3. Histogram of the amount of extinction in magnitudes 
at each position, in bins of 0.05 mag. The blue curve (right) is the 
Aj distribution, green (middle) is the Ah distribution and red 
(left) is the Ak^ distribution. 



The Aks extinction distribution follows approximately 
the same distribution as the Aj and Ah above, in the range 
0.5 ^ Aks ^ 4.5. For Ak^ the low extinction component 
is almost not present, and the higher extinction plateau is 
more pronounced. The main, largest component of the dis- 
tribution is centred on Ak^ =1-1, and the higher extinction 
plateau extends from 1.3 ^ Ak^ ^ 2.5 with two appar- 
ent peaks at Aks = 1.55 and Aks = 2.05. The high ex- 
tinction tail extends from Ak^ = 2.5 out to approximately 
Aks = 4.5. 

Figure 3 shows the three extinction distributions plot- 
ted on the same graph. This shows the differing values of 
the extinction compared to each other, and the similarity 
of the shapes of the three distributions. The three distri- 
butions appear to be a single distribution that is stretched 
at shorter wavelengths {Aj is broader and lower than Aks 
for example). Figure 4 shows the extinction maps of three 
example fields showing the degree of variation within each 
2.5' X 2.5' field. Each pixel in these images corresponds to 
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Figure 4. Maps of the measured extinction in three of the VLT fields in Aj. Each map fills the full extent of the 2.5' X 2.5' VLT 
fields, and each pixel within each map represents the extinction in a region 5" X 5". The colour scale represents variation in the absolute 
extinction in the range: Aj ^ fO. 



one of the grid points separated by 5" as outlined in Section 
4, and the colour scale represents ^ Ah ^ 10. 



5.3 Comparison with non- varying Extinction Law 

In order to test the significance of the variation in the ex- 
tinction law, we re-ran our extinction calculation twice, once 
allowing the value of the extinction law parameter a to vary 
from field to field but not within the field and a second time 
using the mean extinction law parameter for all of the fields, 
a = 2.64 (see Figure 1). To measure the effects, we compared 
how the value of absolute extinction calculated using equa- 
tions 1 changes for the different colour excesses (for example: 
the difference between A,j calculated from {E(J — H)) and 
{E{J — Ks)))- We measured the difference using a xt test 
for each value of A\. 

For the fully varying, V-NICE extinction calculation 
method that we are presenting in this paper, the xt test 
returned values in the range 1.34 — 2.37, indicating that the 
variation between the measures is within one or two a of 
the distribution for almost all cases (Table 1). Allowing the 
extinction law to vary from one field to another, but not 
within the field, we find that the differences in the values of 
absolute extinction calculated from different colour excesses 
result in a xt test returning results of 10.78 — 15.35. This 
illustrates that there is a statistically significant difference 
between the distributions of absolute extinction calculated 
for each wavelength from the two applicable colour excesses 
when keeping the extinction law constant across individual 
fields rather than allowing it to vary within a given field. In 
other words, if the value of a is held constant over the 2.5 
arcmin^ of any one VLT field, different values of the abso- 
lute extinction (i.e.: Aj) will be calculated depending upon 
which colour excess is used (i.e.: {E{J~H)) or {E{J — Ks)}). 
This illustrates that it is important to allow for variation in 
the extinction law on as small a scale as is possible for any 
dataset to derive the most accurate measure of the extinc- 
tion in a given field. 

The final re-run of the extinction calculation used a 
single extinction law, fixed at the average value for all the 
fields {a = 2.64 ± 0.52) as measured in Figure 1, to con- 



vert colour excess to absolute extinction for all 26 of the 
VLT fields. This is similar to the method that has previ- 
ously been used to calculate the extinction throughout the 
Galaxy where variations in extinction law have rarely been 
taken into account. In this case there were obvious differ- 
ences in the calculated values of absolute extinction from 
the different colour excesses, revealed by the xt test which 
returned values of 19.85 — 55.05. 



5.4 Intrinsic effects of filters 

We note that the value of the extinction law parameter a 
that we measure from our data is significantly different from 
the more generally accepted value a ~ 2 (Rieke & Lebofsky, 
1985; Cardelh et al., 1989; Mathis, 1990; Nishiyama et al., 
2006, and others). Here we consider that the method that we 
use to calculate the extinction law parameter may be affect- 
ing the measured value. In equation 3, it is assumed that 
the filters in the telescope system used are infinitely thin, 
whereas in reality they admit light over a range of wave- 
lengths with varying efficiency, the transmission profile. It 
is possible that in sampling a larger bandwidth of the stellar 
spectra than only at the central wavelengths as assumed in 
our methods; the filter transmission profiles are the cause 
of the difference of our measured values of a from previous 
studies. 

To explore this, we simulated the effects of using both 
a narrow idealized filter, as assumed in our extinction cal- 
culation method, and the real transmission functions of the 
VLT-ISAAC fihers as given on the ESQ websiteV To sim- 
ulate the narrow filters we used thin gaussians with peak 
transmission efficiency of 1 and width of 100 A centred on 
the quoted central wavelengths of the VLT-ISAAC filters. 
For the actual filter case, we could only apply the filter trans- 
mission functions as there is no comprehensive data available 
on the ESO website which takes into account the additional 
effects of the atmosphere, telescope and instrument optical 



^ http : //www. eso . org/sci/f acilities/paranal/instruments/ 
isaac/inst/lsaac.img . html 
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Table 1. Results of the tests of the effect of extinction law vari- 
ation on calculated absolute extinction for three scenarios: one 
with a fully variable extinction law (top), one with an extinc- 
tion law fixed for individual fields {middle) and one with a single 

extinction law for all fields (bottom). For each instance, the differ- 
ences in the values of absolute extinction for a single wavelength 
calculated from the two applicable colour excesses arc compared 
with a xi test. The results of these tests are presented against 
the colour excess error used in their calculation. It shows that the 
fully variable extinction law produces the most consistent calcu- 
lation of absolute extinction from different colour excesses and 
that a fixed extinction law produces the worst. 



Colour excess Colour excess 

error used error used x? 



Fully variable extinction law 

Aj{E{J-H)) 2.3675 Aj{E{J-Ks)) 2.2644 

Ah{E(J-H)) 1.3412 Ah {E{H-Ks)) 1.8251 

Aks {E{J-Ks)) 1.7369 A^g (E{H-Ks)) 2.0890 



Siii,;^i(* cxI iiK'l ion law fur twcii [ivUi 

Aj{E(J-H)) 12.0757 Aj{E{.J-Ks)) 10.7796 

Ah{E(J-H)) 12.7166 Aij(E(H-K^)) 15. .3468 

Ak^ {L(J - 12.101)7 11.1556 



tinction law corresponded to variation in the spectral types 
of the stars, earlier type stars resulted in a lower measured 
value of a, and later types resulted in a higher measured 
value of a (within the range given above). 

To average the effects of the different spectral types, 
and to recreate our actual method used to calculate extinc- 
tion, namely sampling a variety of stars over a region of the 
sky, we selected 40-50 random spectral types, adding 5% 
noise to the measured magnitudes in all three bands, and 
calculated an average value for the extinction law for the 
combined colours of the stars. In 1000 repetitions, we mea- 
sured a value of a = 1.81 ± 0.11 (the error being consistent 
with the random noise added to the simulated magnitudes) . 



6 DISCUSSION 

We have identified two significant sources of potential er- 
ror in the calculation of the extinction in the near-infrared. 
First, there may be substantial significant variation in the 
extinction-wavelength power-law relationship parameter a 
(where Ax = C x A^") and in the absolute value of extinc- 
tion Ax, especially in heavily extinguished regions of the 
Galactic Bulge and Plane. Second, the measured values of a 
are influenced by the transmission profile of the instrument 
used to obtain the photometric data on which the extinction 
calculation is based. 



Extinction law a = 2.64 applied to all fields 

Aj{E(J-H)) 19.8532 Aj {E{.J-Ks)) 21.2689 

Ah{E{J-H)) 27.8859 A„ (E{H-Ks)) 55.0459 

Aks {E{J-Ks)) 34.4246 Ak^ {E(H-Ks)) 39.5893 



path on light transmission in VLT-ISAAC. Experience sug- 
gests that these effects will be predominately attenuation of 
the light transmitted through the system and so the filter 
transmission function will have the largest effect on variar 
tion in the measured extinction. We used a sample of stellar 
spectra, G, K and M type giants from the Pickles (1998) 
libraries, and applied extinction to all spectra with a power- 
law wavelength dependence Ax oc A~™. We used a = 2 for 
the extinction applied to the spectra as that is the most 
commonly used value for extinction in the near-infrared and 
it also allows us to compare the results to the values we 
measured with our V-NICE method. We then used the IRAF 
package sjmphot^ to pass the extincted standard star spec- 
tra through the two filter sets described above. We used the 
magnitudes output by the synphot package for the different 
filters to calculate colour excess for the spectral types, and 
from the colour excess calculated the extinction law param- 
eter a using equation 3. 

We found that the case of thin gaussians returns the 
true value of q = 2.00 ± 0.05, as expected. For the actual 
ISAAC filter transmission functions, the values of a mea- 
sured for the different spectral types of the stars were in the 
range 1.75 ^ a ^ 1.85. The variation in the measured ex- 

^ http://www.stsci.edu/resources/software_hardware/ 
stsdas/synphot 



6.1 Extinction calculation 

We have applied our new V-NICE extinction calculation 
method which measures variation on small spatial scales 
in both the extinction law and the absolute extinction to 
fields distributed "randomly" throughout the Nuclear Bulge, 
avoiding known regions of star formation and structures such 
as the Arches and other clusters. As such, these fields of- 
fer a deeper, representative sample of the "general" Nuclear 
Bulge population than it has been possible to obtain from 
previous observations such as 2MASS. As we have demon- 
strated in Gosling et al. (2007), this method produces self- 
consistent measures of the extinction that can recover real 
photometric values for the extinction corrected stars. Un- 
til similar analysis can be undertaken on forthcoming near- 
infrared observations of the entire Nuclear Bulge, using data 
from the UKIDSS and VISTA surveys, this work represents 
the most generalised and representative study of the extinc- 
tion of the overall Nuclear Bulge. 

6.2 Extinction Law 

In all fields we see variation in the measured extinction law 
on scales of ~ 5 — 10", which is also the scale of the gran- 
ularity measured in Gosling et al. (2006). We measure a 
broad range of values of the extinction law parameter, in 
the range ^ a ^5 4. Fitting a gaussian to the extinction 
law parameter distribution, we find that the extinction law 
in the Nuclear Bulge can be described in the most general 
terms of a = 2.64 with a standard deviation of 0.52 (see 
Figure 1). 

We tested the variation in the extinction law parameter 
as measured with VLT-ISAAC, Table 1 shows the effect of 
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the variation on the calculated absolute extinction when us- 
ing different colour excesses. Only in the case of allowing full 
variation of the extinction law are the absolute extinction 
values calculated from different colour excesses the within 
3ct limits. For the cases of a more general extinction law, a 
single value for each field and a single value for all fields, the 
differences introduced by the calculation of absolute extinc- 
tion from different colour excesses are beyond what can be 
explained as statistical variations in the distributions (i.e.: 
differences greater than Scr). Additionally, in Gosling ot al. 
(2007), we demonstrate that the use of a fully variant extinc- 
tion law is the only way to obtain consistent photometry for 
stellar sources in our VLT fields and, by extension, through- 
out the Nuclear Bulge. 

6.3 Absolute Extinction 

Wo sec considerable levels of variation in the measured ab- 
solute extinction on very small scales resulting in overall 
extinction ranges in the near-infrared of 4.5 < Aj ^ 10, 
1.5 ^ Ah ^ 6.5 and 0.5 ^ Aks ^ 4.5. The spread of 
the measured values of absolute extinction are similar for 
all three bands. Each has a small amount of very low lev- 
els of extinction, a large fraction centred on a single value, 
with a "plateau" of increasing values of extinction with a 
final, very small high extinction tail. That these distribu- 
tions are similar, with the shorter wavelength distribution 
"stretched" indicates that all three near-infrared bands are 
measuring the same extincting material distribution. The 
variations we measure are on scales of ~ 5 — 10", the same 
scales at which wc measured variation in the stellar distri- 
bution in Gosling et al. (2006). Wc interpret the distribu- 
tion of different values of extinction as representing the dis- 
tribution of these dense, obscuring clouds identified in our 
previous paper. The largest peak in all three distributions, 
at the lower end of the overall extinction distribution repre- 
sents the "field" value away from these extincting structures, 
the higher extinction "plateaux" represent the extinction in 
the denser regions, the size scale of which we measured in 
Gosling et al. (2006). It is thus important to consider both 
the size of these structures, and their distribution when cal- 
culating the extinction in the Nuclear Bulge as applying a 
"field" value for the extinction could underestimate extinc- 
tion by a factor of 2-3 if the target is close to or in one of 
these dense, complex regions. 

The values of Axg are we measure similar to the values 
that have previously been measured, with the majority of 
the Aks distribution in the region of 1-3; however, due to 
the variable and generally greater value of the extinction law 
parameter o, the values of Aj and Ah in this work are con- 
siderably higher than have previously been measured. Both 
our calculation of the extinction from different colour ex- 
cesses (Section 5.3), and the extinction correction of stars in 
Gosling ot al. (2007) for photometric identification indicate 
that this higher value of the extinction law parameter a and 
the small scale variations are the correct method to use in 
the Nuclear Bulge. The higher value of a means that pre- 
vious works (Catchpole et al., 1990; Schultheis et al., 1999; 
Lombardi & Alvcs, 2001, and others) which quoted values 
of Av based on measures of extinction in the near-infrared 
will have significantly underestimated the values of Ay- As 
an example, the most commonly used conversion factor from 



near-infrared to visual is Av/Aks = 9 (Rieke & Lebofsky, 
1985); more recently Nishiyama et al. (2008) have reported 
that in areas of the outer Bulge, they measured a higher con- 
version factor oi Av /Aks ~ 16. Using the conversion factors 
of Martin & Whittct (1990) and Ricke & Lebofsky (1985) 
to derive the mean relation between Aks and Ay , our work 
would suggest a conversion factor oi Ay /Aks = 28.7l}j y in 
the Nuclear Bulge for a = 2.64 ± 0.52, 3 times higher than 
the most commonly used value. 

6.4 Measuring the Extinction 

We recommend that wherever possible, i.e. if data from at 
least three infrared bands is available for a given location, a 
specific calculation of the extinction be carried out to obtain 
accurate stellar photometry in heavily extinguished regions 
of the sky such as the Nuclear Bulge. Due to the very small 
scale, real variation that we measure in the extinction law 
in this region, it is extremely important that as local a value 
as possible is measured. Using our V-NICE method allows 
for the very small spatial variations in the extinction law pa- 
rameter a as well as the absolute extinction, producing the 
most accurate measure of the extinction at different wave- 
lengths in the near-infrared. 

6.5 Filter transmission effects 

As seen in Section 5.4, the transmission function of the filter- 
set used in observations has an effect on the value of the 
extinction law parameter measured using our method. The 
results suggest that only for an idealized very thin filter pro- 
file will the "real" value of the extinction law be measured 
using this method. For the broader filter profiles of real in- 
struments, the spectral slope of the star being observed will 
influence the value of the extinction law measured. The spa- 
tial variations in the extinction law that wc observe arc not 
affected by this property, only the exact numerical value 
measured for a. The result of this is that the extinction 
law parameter a and the absolute extinction measured from 
a sot of observations using a unique telescope-instrument 
setup can only be applied to data from other observations 
using the same telescope-instrument combination without 
the need to account for the filter transmission. If the extinc- 
tion correction values measured from one dataset are applied 
to data from observations obtained using other telescope- 
instrument combinations, an error will be induced by the 
differing filter transmissions. 



7 CONCLUSIONS - A NEW EXTINCTION 
MEASURE 

Wc have studied the extinction in the near-infrared in a sam- 
ple of representative regions of the Nuclear Bulge based on 
observations with VLT-ISAAC. We have found a high degree 
of spatial variability in both the extinction law parameter a 
as well as the absolute extinction in these regions. We have 
used a now method, a variable near-infrared colour excess, 
V-NICE, to measure these variations. Only when the value 
of the extinction law parameter a as well as the absolute 
extinction is allowed to vary can self-consistent values for 
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the extinction at different wavelengtlis in the near-infrarod 
be calculated. 

For our observations, using VLT-ISAAC, the peak of 
the extinction law parameter a's distribution for the Nu- 
clear Bulge is at a higher value than previously thought, 
described by a gaussian with = 2.64 and a = 0.52. We al- 
lowed the extinction law to vary when calculating the abso- 
lute extinction at different wavelengths from colour excesses 
and show that only when this variation is allowed can self- 
consistent levels of extinction be recovered. The variations of 
both the extinction law parameter a and the absolute extinc- 
tion are seen on scales as small as 5" . We measure extinction 
in the near-infrarod of 4.5 ^ Aj ^ 10, 1.5 ^ Ah ^ 6.5 and 
0.5 ^ Aks ^ 4.5. The higher overall value of the extinc- 
tion law parameter a results in a higher conversion factor 
between Aj, Ah and Aks- This suggests that previously 
derived levels of Ax in the Nuclear Bulge from near-infrared 
extinction measures may need to be reconsidered. 

We have also found that the extinction law parameter 
a has a dependency on the transmission function of the in- 
strument. We simulated the effects of the instrument filter 
transmission profiles on the measured value of the extinc- 
tion law parameter a and showed that in the case of ide- 
alized narrow filters, our method recovers the true value of 
the extinction law, but for the real ISAAC filter transmis- 
sion profiles, the result measured is ~ 10% below the true 
value. 

The effects of the small scale variation in the extinc- 
tion law as well as the absolute extinction axe the domi- 
nant source of errors when correcting for extinction in the 
Nuclear Bulge. Furthermore, wc suggest that the potential 
for small-scale physical variations in the extinction should 
be considered when performing photometry in any signifi- 
cantly extinguished area of the Galactic Bulge or Plane. The 
method presented here is sufficient to obtain accurate stel- 
lar photometry for multi-colour data from any near-infrared 
telescope-instrument combination, and can be used to cor- 
rect for extinction to allow identification of stellar types by 
comparison to photometric standards. The exact extinction 
values measured will be relative to the telescope-instrument 
combination. Unless the absolute value of extinction law is 
required, or the measured extinction will be used to cor- 
rect the colours of stars from datasets from other combinar- 
tions of telescope and instrument, the effects of the filter 
transmission effects discussed here need not be considered. 
However, if the measured extinction law value is used to 
reddening-correct data from different telescope-instrument 
combinations, the user should be aware of the errors that 
will be induced by the different filter transmission profiles. 

Due to these two effects, we suggest that the use of a 
generic "universal" IR extinction law is not appropriate in 
many cases, especially when observing regions of the Galaxy 
known to suffer significant extinction. We recommend that 
whenever possible, a local extinction calculation should be 
undertaken as this provides the most accurate possible pho- 
tometry for Bulge stars (Gosling et al., 2007). We suggest 
that the measure of granularity from Gosling et al. (2006) be 
used to determine the scale at which the extinction is vari- 
able and this be used to determine a limit to the area over 
which stars can be sampled in calculating the extinction. 
The methods we present here are an effective and accurate 
way of calculating the extinction for a dataset to correct 



the extinction towards sources in the Nuclear Bulge in the 
near-infrared. 
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